arXiv:1509.07176v2 [quant-ph] 13 Dec 2015 


Position-Momentum Bell-Nonlocality with Entangled Photon Pairs 


James Schneeloch/’^ Samuel H. Knarr/’^ Daniel J. Lum/’^ and John C. Howelh’^ 

^Department of Physics and Astronomy, University of Rochester, Rochester, NY 14627 
^ Center for Coherence and Quantum Optics, University of Rochester, Rochester, New York 14627, USA 
^Air Force Research Laboratory, Information Directorate, Rome, NY 1344^, USA 

(Dated: December 15, 2015) 

Witnessing continuous-variable Bell nonlocality is a challenging endeavor, but Bell himself showed 
how one might demonstrate this nonlocality. Though Bell nearly showed a violation using the CHSH 
inequality with sign-binned position-momentum statistics of entangled pairs of particles measured 
at different times, his demonstration is subject to approximations not realizable in a laboratory 
setting. Moreover, he doesn’t give a quantitative estimation of the maximum achievable violation 
for the wavefunction he considers. In this article, we show how his strategy can be reimagined using 
the transverse positions and momenta of entangled photon pairs measured at different propagation 
distances, and we find that the maximum achievable violation for the state he considers is actually 
very small relative to the upper limit of 2^/2. Although Bell’s wavefunction does not produce a 
large violation of the CHSH inequality, other states may yet do so. 


I. INTRODUCTION 

Though there are many predictions of quantum me¬ 
chanics that disagree with classical intuition, perhaps the 
most striking finding is that quantum mechanics predicts 
the violation of Bell inequalities, ruling out the possibil¬ 
ity that correlations between distant events can always 
be explained by shared information in the past. Violating 
Bell inequalities is more than experimental metaphysics; 
Bell-nonlocal entangled systems (i.e., those whose statis¬ 
tics cannot be described by a Local Hidden Variable 
(LHV) model) can be used to prove secure key rates 
in device-independent quantum key distribution [T]. Of 
particular significance in quantum information protocols 
is the demonstration of Bell nonlocality in continuous- 
variable (CV) systems; their high dimensionality offers 
the possibility of transmitting much more information 
with individual particles than what their spins or polar¬ 
izations can convey. 

To demonstrate Bell-nonlocality [5] in continuous ob¬ 
servables, one must violate a continuous-variable Bell in¬ 
equality. By some standards, a fully general (i.e., contex¬ 
tual) continuous-variable Bell inequality does not yet ex¬ 
ist, since we cannot deduce Bell nonlocality (yet) simply 
from knowing the first and higher-order moments of con¬ 
tinuous observables for different measurement settings. 
However, there has been much research into demonstrat¬ 
ing continuous-variable nonlocality by examining low¬ 
dimensional observables derived from the continuous ob¬ 
servables of interest (examples include pseudo-spin ob¬ 
servables [3] and parity observables [4]). In addition, 
there have been efforts in demonstrating continuous- 
variable Bell nonlocality with the statistics of discrete 
functions (e.g., binning functions) of continuous observ¬ 
ables [5Hin! (in these cases, field quadratures). Indeed, 
since an LHV model for a continuous-variable joint prob¬ 
ability distribution is also a model for its consequent 
statistics, violating the CHSH inequality with such (e.g., 
discrete) statistics will demonstrate the nonlocality of the 
underlying continuous observables. 


In fact, one of the first historical attempts at showing 
that continuous-variable nonlocality was possible was de¬ 
veloped by Bell himself mm, who showed that with 
sign-binning, there were wavefunctions that may violate 
the Clauser-Horne-Shimony-Holt (CHSH) inequality |13j . 
Ironically, he also showed that the position-momentum 
statistics of the maximally entangled Einstein-Podolsky- 
Rosen (EPR) state cannot violate a Bell inequality, since 
they admit an explicit LHV model. 


In spite of theoretical demonstrations that it is pos¬ 
sible to demonstrate CV Bell nonlocality with sign bin¬ 
ning, there have been no experiments to date that give a 
successful demonstration of Bell nonlocality for continu¬ 
ous observables (though there has been progress in doing 
so for field quadratures [MH])- In this article, we show 
how one might plausibly demonstrate CV Bell nonlocal¬ 
ity with the transverse spatial statistics of entangled pho¬ 
ton pairs. Euthermore, we show how the nonlocal state 
Bell considers actually gives a very small violation of the 
CHSH inequality. 


To those familiar with the transverse spatial statistics 
of highly entangled photon pairs, it may seem impossi¬ 
ble to demonstrate nonlocality in that degree of freedom. 
Indeed, it is reasonably popular to approximate the joint 
transverse spatial amplitude of such photon pairs as a 
Double-Gaussian [niin] function. Such wavefunctions, 
have Gaussian Wigner functions, which are non negative, 
and so admit an LHV model (as shown later). However, 
there are multiple states of entangled photon pairs (see 
Fig. 1) whose Wigner functions have significant regions 
of negative values. Though this does not, by itself guar¬ 
antee nonlocality, it provides sufficient motivation to see 
if there are states other than the one Bell considers that 
might give a more substantial violation. 
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A. Continuous-Variable CHSH Inequality with 
Sign Binning 


pendent of conditioning. With a little extra algebra (see 
[13)1. one can show that the inequality: 


In a narrow sense, all Bell inequalities descend in one 
form or another from the same model of contextual [16j 
Local Hidden Variables (LHVs) that gave rise to the orig¬ 
inal CHSH inequality |13| . As such, many current tests 
of Bell nonlocality rely on one form or another of the 
CHSH inequality [T7] . 

To begin, an LHV model for, say, the measurement 
outcomes of the positions xi and X 2 of a pair of parti¬ 
cles is one where the correlations between Xi and X 2 are 
explained (if not also determined) by a complete knowl¬ 
edge of every piece of information that could possibly be 
conveyed to each particle (at or below the speed of light) 
from points in the past. To describe these pieces of in¬ 
formation in a general way, we assign the variable(s) A. 
These models are local in that A (at the very least) en¬ 
codes all information in the intersection of the past light 
cones of both particles. If the position correlations can be 
explained by all local information A, then the joint prob¬ 
ability density p(xi,X 2 ) can be expressed in the following 
form: 

p{xi,X 2 ) = j dX {p{X)p{xi\X)p{x 2 \X)). (I) 

The CHSH and many other Bell inequalities, are math¬ 
ematical consequences of LHV models of this form. By 
violating the CHSH inequality, we rule out the possibil¬ 
ity that the joint probability density can be expressed 
in this form. If we can show that the joint probability 
density is not expressible in this way, it then follows that 
knowing every piece of information in the past light cones 
of both particles cannot explain the correlations between 
them. This opens up new possibilities where there is ei¬ 
ther information outside both particles’ light cones that 
will explain the correlations (implying a nonlocal uni¬ 
verse), or that there is no information at all that will 
explain these correlations (implying a non-deterministic 
universe). As quantum metaphysics, these are interesting 
questions, but well beyond the scope of this work. 

The derivation of the CHSH Bell inequality is based 
both on the essential LHV model, and on the measure¬ 
ments considered having a bounded (though possibly 
continuous) spectrum of possible results (TB]. To keep 
things general, let a; be a real-valued random variable, 
and f{x) be a function whose range is in the closed set 
[—1,1]. Let a, and a' be two possible measurement set¬ 
tings for one party, and let (3 and /?' be two possible 
measurement settings for a second party. If the position 
statistics can be described by an LHV model Q, then 
the correlation {f{xi)f{x 2 ))a ,0 must be expressible as 
the form: 

{f{xi)f{X2))a,0 = J dX p{X){f{xi))a,x{f{x2))p,\, ( 2 ) 

for some A. Next, since for all x,f{x) G [—1,1], all ex¬ 
pectations (fix)) must fall within the same range, inde- 


\{f{xi)f{x 2 ))c,fi - {f{xi)f{x 2 ))a,fi'\ 

T {{f{xi)f{x 2 ))a',P + {f{xi)f{x 2 ))a',fi') < 2, (3) 

must be valid if an LHV model of the form Q exists. 
Indeed, when x has a binary spectrum (i.e., {—1,-|-1}), 
we see how the CHSH inequality as it applies to polariza¬ 
tion measurements [MS] is a special case of this more 
general formulation. 


B. Bell’s Wavefunction and CHSH Violation 


What makes Bell nonlocality in continuous variables 
especially difficult is that many useful approximations 
to continuous-variable quantum states have positive- 
definite Wigner functions, and so explicitly admit an 
LHV model. To see how this works, we consider the 
Wigner function W{xi,X 2 ,ki,k 2 ) for a pair of particles 
(in one spatial dimension) with position observables Xi 
and X 2 , and momentum observables ki and /c 2 - The joint 
probability density p{xi,X 2 ) can be expressed in terms 
of the Wigner function: 


p{Xi,X2) 


dkidk 2 W{xi,X 2 ,ki,k 2 ), (4) 


where: 

Wixi,X2,kuk2) = ^ Jj 

X ^(^1 + - |^,fc 2 - y). ( 5 ) 

Note that if the correlations admit an LHV model, then 
p{xi,X 2 ) has the form seen in equation Q. Because of 
this, when the Wigner function happens to be positive- 
definite (obeying the form of a probability density), we 
can immediately describe p{xi,X 2 ) with an LHV model; 
the hidden variables A would be the arguments of the 
Wigner function, while the probability densities in paren¬ 
theses would constitute the Wigner function itself. 

Since wavefunctions with non-negative Wigner func¬ 
tions admit an LHV model, we know that such states 
will never violate Bell inequalities based on those mea¬ 
surement statistics. In spite of these difficulties, there 
are relatively simple wavefunctions that can violate a 
continuous-variable Bell inequality. Indeed, Bell pro¬ 
vided a specific example in |12] : 

_ (3:i+3:2)^ _ 

^^^^\xi,X 2 )=N{{xi-X 2 f-Sal)e ""+ e , 

( 6 ) 

where iV is a normalization constant. This “Bell- 
Violating” wavefunction is especially convenient because 
it is separable in terms of the rotated (orthogonal) coor¬ 
dinates: 


Xi -X- X2 

^7^ 


Xi — X2 

~7pr 


x+ = 


X- = 


(7) 
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FIG. 1: (Color online) On the left (a) is a plot of the Wigner function W_ {x-,k-) obtained from direct 
calculations of the biphoton amplitude from Spontaneous Parametric Down-Conversion (SPDC), a popular source of 
entangled photon pairs (as shown in [I1I11I23])- We note that although the biphoton Wigner function can be 
approximated as a Gaussian, there are significant regions of negativity. On the right (b) is a plot of Bell’s Wigner 
function k-) (111 for ct_ chosen to match position moments of the SPDC biphoton state. The values 

below the level (thick red) contours are negative. 


As a result, the Wigner function for is separable 

this way as well, i.e., 

W^^^'>{xi,X 2 ,kuk 2 ) = {x +, /c+) {x- fc_), 

( 8 ) 

where 

Wix+,k+)^^ II + 

(9) 

and W{x-, k-) is similarly defined for x- and fc_. Now, 
W'^^\x+,k+) is a Gaussian: 

kl 

1 -7-p- 

Wi^^\x+,k+) = -e , ( 10 ) 

TT 

but since ip^^'^'>{xi,X 2 ) has both quadratic and Gaus¬ 
sian factors depending on x-, W^^\x-,k-) is a more 
elaborate function with negative values: 

W[^'^\x_,k_) = — ^(xt+2xW_{-5 + ik^_a^_) 
ll7r(7_ \ 

-fcri(ll-f Sfc^cr^ -P 16/clcri))e”^e . (11) 

In Fig. lb, we plotted lFl^^^(a;_, /c_) to show where the 
negative values occur. 

To demonstrate Bell nonlocality for continuous vari¬ 
ables (in this case, the positions xi and X 2 of a pair of 
particles), Bell examined (for different measurement set¬ 
tings) the signs of the position measurement outcomes 
(-1-1 for X > 0, and —1 for x < 0). He then took 
for his correlation measure the mean of the product of 
these signs. Since these mean products are bounded be¬ 
tween — 1 and 1, these statistics can be readily used in 


the CHSH inequality to test for nonlocality in continu¬ 
ous variables (e.g., /(x) in could be Bell’s piecewise 
sign function). To examine these correlations for differ¬ 
ent measurement settings, he treated the pair of systems 
described by the wavefunction as though they were 
free non-interacting particles, and time evolved their joint 
state to different points in time ti and ^2 according to 
the free particle Schrodinger equation. Calculating the 
sign correlation at different pairs of times allowed him to 
show by example how one might violate the CHSH Bell 
inequality with continuous-variable states. 


II. ADAPTING BELL’S APPROACH TO 
MEASUREMENTS OF PHOTON PAIRS 

Although Bell showed theoretical violation of the 
CHSH inequality by propagating his wavefunction with 
the free particle Hamiltonian to different times, one may 
violate the CHSH inequality experimentally by consid¬ 
ering Bell’s wavefunction as describing a biphoton trans¬ 
verse position amplitude (say, for signal and idler photons 
in degenerate spontaneous parametric down-conversion 
(SPDC)), and considering the sign correlation measure¬ 
ment performed at different propagation distances, as 
opposed to different times. The justification for this 
comes from the fact that the paraxial Helmholtz equa¬ 
tion (governing the propagation of the individual pho¬ 
tons) is mathematically identical (variable names aside) 
to the Schrodinger equation for a free particle moving in 
two transverse dimensions, e.g., 

d^A d'^A dA ,2m d'i’ 

9x^ dy'^ * ^ dz 9x^ dy"^ * dt' 

( 12 ) 
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3 ) CHSH Violation vs. Propagation Distances 
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b) CHSH Violation VS. Misalignment 
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FIG. 2: (Color online) On the left (a) is a contour plot of the violation of the CHSH inequality for different 

propagation distances zi and Z 2 - These plots were calculated using Bell’s wavefunction with = 1mm, 

cr_ = 0.01mm, and a wavelength of 650nm for each photon. On the right (b) is a contour plot of the violation of the 

CHSH inequality at the optimal values of Zi and Z 2 (—3 and +8 mm respectively) for different misalignments of the 

detectors, dxp = dxi + dx 2 is a parallel position misalignment (both shifted the same direction), while 

dxm = dxi — dx 2 is an anti-parallel misalignment (both shifted in opposite directions). A negative value indicates a 

violation of the CHSH inequality. 


where A(x, y, z) gives the spatial dependence of the am¬ 
plitude of either the signal or idler electric field; and kp 
is the wavenumber of the pump electric field (i.e, twice 
the signal or idler wavenumbers). 

Using his Wigner function W^^^'>{xi,ki-,X 2 ,k 2 ), in the 
approximation that cr+ ^ cr_ (so that W^'^\x+,k^) 
may be neglected as a constant factor). Bell showed a 
qualitative violation of the CHSH inequality (§ for op¬ 
timum measurement settings. However, if we are to use 
this relationship (12) to get a quantitative violation, we 
need to show that a violation is possible without such 
approximation, as the diffraction of light implies this ap¬ 
proximation only works for a narrow range of distances. 
Indeed, in Bell’s own treatment, the spreading of the en¬ 
tangled wavepacket over time implies his approximation 
(and subsequent result) is only valid for a narrow range 
of time settings [21] • 


A. Quantitative Violation of CHSH Inequality 
with Bell’s Wavefunction as a Biphoton Amplitude 


values of tT+ = 1mm and ct_ = 0.01mm, we were able 
to show that we could violate the CHSH inequality 
but only by at most a small amount (i.e., our maximum 
value was 2.00041, though with a numerical error bound 
less than 10“®) (see Fig. 2). In order to realize such 
a minute violation, we would need approximately 10® 
coincidence counts to barely resolve each of the sixteen 
probabilities in the CHSH inequality to an uncertainty 
of ±1.0 X 10“^, and that is assuming any systematic er¬ 
ror due to misalignment (see Fig. 2b) or other factors is 
insignificant. We also performed similar calculations for 
the biphoton state created from type-1 SPDC, but any 
violation was inconclusive as the uncertainty in the re¬ 
sults from these highly oscillatory probability densities 
was larger than the obtained violations. Although this 
underscores the difficulty in a successful demonstration of 
position-momentum Bell-nonlocality, it is still useful to 
consider what sort of experimental setup might be used 
in a successful demonstration once a more suitable state 
is discovered. 


In order to see whether it is empirically feasible to 
violate the CHSH inequality for position and momen¬ 
tum with sign binning, we need to be able to violate this 
inequality without taking such limits (e.g., a+ ^ ct_). 
By performing numerical calculations of the Fresnel- 
propagated biphoton field, we find that with a biphoton 
field resembling Bell’s wavefunction, and using typical 


B. Imagining an Experimental Demonstration with 
Photon Pairs 

The idealized setup we consider is as seen in Fig. 3. 
A pump laser and nonlinear crystal serve as a source of 
entangled photon pairs though other sources may be sub¬ 
stituted. The pump light is filtered out, and the photon 
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FIG. 3: (Color online) Here is a diagram of an idealized 
setup to violate the CHSH inequality for position and 
momentum with sign-binning. A laser and a nonlinear 
crystal (NLC) serve as a source of entangled photon 
pairs, whose state is determined by both the particulars 
of the laser and the geometry of the NLC. We would 
then image the exit face of this photon-pair source onto 
pairs of avalanche photo-diodes (APDs) with a 4F 
imaging system. By translating the APDs along the 
optic axis, we may use coincidence counts to measure 
the transverse sign correlations as a function of a 
variable propagation distance in each arm. 


pairs are separated by a 50/50 beamsplitter into separate 
arms. Each arm contains a 4F imaging system, where 
spatially resolving pairs of photon detectors are placed 
in the image planes conjugate to the exit face of the pho¬ 
ton pair source. In addition, these detectors are placed on 
translation stages, allowing us to measure the sign corre¬ 
lations for different propagation distances in the different 
arms. 

The reason we would need a 4F imaging system in 
each arm is that the biphoton field formed in those con¬ 
jugate planes would be identical to the field just as it exits 
the source (though reflected and subject to the paraxial 
approximation). Without it, negative propagation dis¬ 
tances would be impossible to realize, as there are no 
photon pairs before the source [25]. 

Though the 50/50 beamsplitter halves the collection 
efficiency of photon pairs, it does not alter the trans¬ 
verse spatial statistics of the photon pairs. Alternatively, 
if we use a source of entangled photon pairs of orthog¬ 
onal polarizations (as with type-2 SPDC), a polarizing 
beamsplitter can better separate the pairs, allowing for a 
larger collection efficiency. Here, we are interested only in 
an idealized setup that may be improved upon in future 
developments. 


C. Engineering a Suitable Biphoton Amplitude 

Once a more suitable biphoton amplitude is discov¬ 
ered, it is important to consider how such a state can 
be created experimentally. Here we consider engineering 
a biphoton state from a source of SPDC photon pairs. 
The biphoton amplitude in SPDC is determined by two 
major factors. The first factor is the pump spatial pro¬ 
file, which we may control with a spatial light modula¬ 
tor, and standard optical components. The second fac¬ 
tor is how the second-order nonlinear coefficient varies 
over the length of the nonlinear crystal (what we call the 
longitudinal nonlinearity profile). For standard nonlin¬ 
ear crystals, this nonlinearity profile is a constant top- 
hat function over the propagation distance z (being one 
value within the crystal, and zero outside the crystal). 
As discussed in [2B], the transverse biphoton amplitude 
(in momentum space) is related to the Fourier transform 
of the nonlinearity profile of the crystal. If we could con¬ 
tinuously vary the nonlinearity within the nonlinear crys¬ 
tal, while controlling the pump spatial profile, we could 
exactly reproduce Bell’s wavefunction among many oth¬ 
ers as a transverse spatial amplitude for photon pairs in 
SPDC. 

Though we cannot continuously vary the nonlinearity 
profile of the crystal, it is possible (as also shown in [26|) 
to use a periodically poled nonlinear crystal, and adjust 
the duty cycle (i.e., the fraction of positive to negative 
poling within each poling period) as a function of the 
crystal length to get a biphoton amplitude closely resem¬ 
bling the state we want. The high fidelity required for the 
approximating biphoton wavefunction in the periodically 
poled nonlinear crystal make this approach challenging, 
though not outside the realm of possibility. 


III. CONCLUSION 

Historically, the CHSH inequality and other Bell in¬ 
equalities have largely been used for finite-dimensional 
discrete observables. However, Bell himself showed how 
one might use the CHSH inequality to witness the non¬ 
locality of continuous observables. Since a local hidden 
variable model for a continuous-variable joint probability 
distribution implies a similar model for any statistics de¬ 
rived from that distribution, violating the CHSH inequal¬ 
ity with sign binned statistics demonstrates the nonlocal¬ 
ity of the underlying continuous observables. Although 
Bell showed a qualitative demonstration of position- 
momentum Bell-nonlocality, we showed that the state 
he considers gives a small quantitative violation, though 
other states may yet do much better. In addition, we 
showed how his approach can be re-envisioned as an ex¬ 
periment measuring similar correlations between entan¬ 
gled photon pairs, and discussed important challenges 
and issues, which when overcome, lead to a successful 
demonstration of position-momentum Bell-nonlocality. 
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